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The baryon-baryon interactions with strangeness S = —2 with the flavor SU(3) breaking 
are calculated for the first time by using the HAL QCD method extended to coupled 
channel system in lattice QCD. The potential matrices are extracted from the Nambu- 
Bethe-Salpeter wave functions obtained by the 2 + 1 flavor gauge configurations of CP- 
PACS/JLQCD Collaborations with a physical volume of (1.93 fm) 3 and with m^/mx = 
0.96, 0.90, 0.86. The spatial structure and the quark mass dependence of the potential 
matrix in the baryon basis and in the SU(3) basis are investigated. 

Subject Index B38, B64, D32, D34 


1. Introduction 

Studying the baryon-baryon ( BB ) interactions in strangeness S = —2 channel is an impor¬ 
tant step to understand hypernuclei such as double-A hypernuclei and 5 hypernuclei (see 
e.g. [1-3]) as well as exotic hadrons such as the iT-dibaryon [4, 5]. Moreover, the hyperon 
superfluidity in the core of the neutron stars is intimately related to the hyperon-hyperon 
interaction in the S = —2 channel [6]. Due to the limited experimental data, however, the 
BB interactions in the S = —2 channel are far from being realistic even under the constraints 
from the approximate flavor SU(3) symmetry. In addition, there are nearby two-baryon states 
in the S = —2 channel (e.g. AA and Nr,) so that the coupled-channel treatment is essential 
for studying the S = —2 system. 

Recently, the BB interactions in the flavor SU(3) limit have been studied systematically 
in full QCD simulations on the lattice by the HAL QCD method (reviewed in [7]) for several 
different masses of the pseudo-scalar meson mps = 470 ~ 1170 MeV (see [8] and references 
therein). In this case, all the two-baryon thresholds are degenerate and the classification 
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of the BB interactions in the flavor basis is applicable. An extension of the HAL QCD 
method beyond the inelastic threshold, which is relevant for the BB interactions with the 
flavor SU(3) breaking, has been also proposed to treat coupled channel systems[9, 10]. The 
main purpose of this paper is to study the BB interactions in the S = — 2 channel with the 
explicit SU(3) breaking on the basis of the coupled channel HAL QCD method developed 
in our previous works [9, 10]. 

We note here that the hyperon-nucleon scattering length away from the SU(3) symmetric 
limit was first evaluated by quench QCD simulation in [11] and later by full QCD simu¬ 
lations in [12, 13], where the Liischer’s finite volume method was utilized. On the other 
hand, the hyperon-nucleon potentials, which provide much more information than the scat¬ 
tering lengths, were derived through the equal-time NBS wave functions with the HAL QCD 
method in [14-16]. The present study can be regarded as a coupled channel generalization 
of our previous works below the inelastic threshold. 

This paper is organized as follows. In Sect. 2, we review the coupled channel approach to the 
the BB interactions by the HAL QCD method in lattice QCD. In Sect. 3, we define baryon 
operators and the baryon states. In Sect. 4, numerical setup on the lattice is summarized. 
In Sect. 5, we present our numerical results of the BB potentials. Summary and conclusions 
are given in Sect. 6 . 

2. Coupled channel BB potentials 

In this section, we briefly review the coupled channel HAL QCD method [9, 10] applicable 
to the inelastic scattering that a\ + —>• b\ + 62 , where (ai, 02 ) Q ( 61 ,Q). 


2.1. Formalism 

We first define equal-time Nambu-Bethe-Salpeter (NBS) wave functions with the total energy 
Wi as 


Wit 

Tp b wS^ e ~ Wit 


1 



VZbi \J Zb 2 


£*< 0 


£*<o 


B ai (x + r,t)B a2 (x,t)\B = 2, Wi), 
B bl (x + r,t)B b2 (x,t)\B = 2,Wj), 


( 1 ) 


where B Cj (x,t) with c = a, b and j = 1,2 denotes a local composite operator for baryon 
B c . with its wave-function renormalization factor yJZ Cj . The state |B = 2, HQ stands for a 
QCD asymptotic in-state with baryon number 2 and energy W\. In the present exploratory 
studies, we assume V ZalV Za2 — V ZblV Zb2 which implies that the flavor SU(3) breaking in 
the wave-function renormalization factor is not sizable in the present setup. The validity of 
this assumption is left for future studies. 

In the asymptotic region at long distance, these NBS wave functions satisfy free 
Schrodinger-type equations as 

+ b) = 0> < 2 > 


where the corresponding asymptotic momentum k\ in the center-of-mass (CM) frame is 
defined through the relation, 


Wi = 


\j™\ + iK ) 2 + 


sJm 2 C2 + (Af) 2 , 


( 3 ) 
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with m Cj being the mass of the baryon B c , and the reduced mass p c is given by 1/ p c = 
1 /m Cl + l/m C2 . On the other hand, in the interaction region at short distance, we have 

K c (r, Wi) = (^0- + r w Jr) + 0, (4) 

from which we define the energy-independent non-local potential matrix as 



This is an extension of the HAL QCD definition for the potential to the coupled channel 
case [9]. To handle the non-locality of the potential, we introduce the derivative expan¬ 
sion as U(r, f 1 ') = (14,0(r) + 4nlo(a) + ■ ■ ■ )<5(r — 4), where N n LO term is of 0(V n ). At low 
energies, a good convergence of derivative expansion has been confirmed for the NN case [17]. 

2.2. Extraction of potential matrix 

In the leading order of the derivative expansion of the non-local potential, eqs. (4) and (5) 
can be written as a coupled channel form of the Schrodinger equation for two independent 
channels a and b, 

( (Ef-HS)r Wt (r) V\(r) V\(r) \ ( ^(r) \ 

^ (E b — Hq) ifyy. (f) ) \ V b a {r) V b b {r) ) ^ if b m {r) ) 1 J 

where the kinetic energy and the free Hamiltonian for channels c = a,b are given by Ef = 
and Hq c = — respectively. 

Two pairs of NBS wave functions, , tfyy. }i=i,2 are necessary to extract the local poten¬ 
tial matrix from the above coupled channel equation. In the infinite volume, we can have 
two states, |a, W) and | b, W) with a given energy W, which are connected to the asymptotic 
scattering states if W is larger than m ai + m a2 and m bl + m b2 . This implies that two nearby 
eigenstates, |B = 2, W\) and |B = 2, Wf) with W\ — W 2 = 0(L~ 2 ), exist even for finite vol¬ 
ume. Suppose that W\ < W 2 are two lowest energies of two baryons in the finite volume. 
By using the wall-source operators Z a (t) = {B a2 B ai )(t) andlfc(t) = (B b2 B bl )(t) 1 , the states 
|.B = 2, W\) and \B = 2, Wf) are created as 

Z c (0)|0) = C cl \B = 2,W!) + C c2 \B = 2,W 2 ) + ■ ■ ■ , (7) 

where the coefficient matrix C c j can be determined from two-baryon correlation functions. 
We then define optimized source operators as 

O «)-(£■ f Y(fit)’ (8 > 

y z W 2 (t) j y c b \ c b i j y i b [t) j 

so that four point (4-pt) function Ff w (r, t ) at large t behaves as 

F x Wi (r,t) = {0\B Cl (x + r,t)B C2 (x,t)l Wi (0)\0) + 0 (e~ W3t ) (9) 

for i = 1,2 and c = a, b, where W3 corresponds to the 3rd state satisfying W\ < W 2 < W3 < 
Wj> 4. By using these 4-pt functions, the coupled channel potential matrix can be determined 


1 A detailed definition of the wall source operators will be given in Sect. 4. 
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for sufficiently large t, where the states with W{> W\, W 2 ) can be neglected in the above 
4-pt functions. As the volume increases, however, the spectrum becomes denser and two 
lowlying states W\ and W 2 cannot be isolated unless extremely large t is achieved. This is 
why we need an improved method in practice, as explained in the next subsection. 


2.3. Time-dependent method 

The improved method to extract the potentials without using the ground state saturation 
has been proposed in Ref. [18] in the case of the single channel. In this subsection, we extend 
this method to the coupled channel case. 

We first introduce the normalized 4-pt correlation function R defined as 




FfJrJ) 


exp [—(m Cl +m C2 )t\ 


^ ( 


re 


-A Wft a W; 




( 11 ) 


where A W? = Wj — m Cl — m C2 and A ^ = (W r j|X ( i(0)|0). The 4-pt function F^ d (f,t ) here is 
defined through the original wall-source operator 2^(0) instead of 1^(0). The ellipses in 
Eq.(ll) denote inelastic contributions from channels other than a and b. 

In the non-relativistic approximation valid at low energies, AIT? ~ Ej, we can replace the 
kinetic energy term in the equation with the time derivative as 

- (r - i) ~ Y, Ejtw, (Oe- AWl, Aj ’, (12) 

3 

with which we obtain the Schrodinger type equation, 


5 H 


R c Xd (r,t) = d 3 r'U c e (ry)A c e R e Id (7,t), 


(13) 


where A c e = exp[— (m ei + m e2 )t\/ exp[— (m Cl + m C2 )t}. If we go beyond the non-relativistic 
approximation, higher-order time-derivatives appear, which we will not consider in this 
paper. Expanding U in terms of derivatives again, the leading order coupled channel 
potentials can be obtained as 


f V\{r) V\(r)A\ \ 

^ V b a (r)A b a V b b (r ) J 


(-Wt-H$)R a Xb (Et)\ ( RliEt) RZ b (r,t) V 

V (-l-^o)^(^*) (-§i- H o)lA h m) \ R b Ia (rA) R^(f,t) ) ' J 

Extension of this formula to three channels is straightforward. For eq.(14) to work, two 
independent source operators T a and I b are needed, while no optimization is required. Note 
that isolation of each eigenstates is not necessary in this method [18]. Only the constraint 
is to keep moderately large t so that other channels having larger threshold energies than a 
and b can be suppressed. In the following, we employ this improved method in our numerical 
calculations. 


4/17 





Table 1 

Summary of channels with S = —2. 



Channel 

Baryon-pairs 

SU(3) multiplets 

1 = 0 

'So 

AA, (NE) o, (SS)o 

1, 8 a , 27 


3 5i - 3 D 1 

(NZ) o 

8 a 

I = 1 

1 So 

(NS),, (AS), 

8 S , 27 


to 

Cp 

1 

to 

b 

(NS) 1; AE, (EE)i 

8a, 10, 10 

1 = 2 

'5 0 

(SS) 2 

27 


3. Strangeness S = — 2 two-baryon system 

We emply the following interpolating operator for octet baryons, 

B a (x) = e abc {ql(x)C^q b (x))q ca (x) (15) 


with the Dirac index a, which represents the spin of the octet baryons. Denoting quark- 
flavors as q = u,d,s for ”up”, ’’down” and ’’strange”, respectively, the flavor structures of 
baryons are given in terms of the isospin multiplets as 


5 = 0 I = 1/2 : p = [ud]u , n = [ud]d 

5 = —1 1 = 1 : E+ = — [us]u , E° = — ([ds]u + [us\d)/y/2 , E - = — [ds\d 

S = — 1 1 = 0 : A. = {[sd}u + [us\d — 2[du\s) / 

S = — 2 I = 1/2 : 5° = [sti]s , S - = [sd]s 


Considering the Fermi-Dirac statistics of two baryons, the allowed combinations for the 
S = — 2 system are given in Table 1, where I z = 0 components are given as 


(EE), 


(NZ)i 


J\ (E+E- + S"E+ - E°E°) , 1 = 0 

(E+E- - S-E+), 1 = 1 

J\ (E+E” + E"E+ + 2E°E°) , 1 = 2 

(pE~ — nS°) , 1 = 0 

\ \[\ + nS °) 5 1=1 


(17) 


(18) 


4. Numerical simulations 

We employ 2 + 1-flavor full QCD gauge configurations from Japan Lattice Data 
Grid(JLDG)/International Lattice Data Grid(ILDG) [19]. They are generated by the CP- 
PACS and JLQCD Collaborations [20] with the renormalization-group improved gauge 
action and the non-perturbatively 0(a ) improved Wilson quark action at j3 = 6 /g 2 = 1.83 
(corresponding lattice spacing in the physical unit, a = 0.1209 fm [21]) on a L 3 x T = 
16 3 x 32 lattice (corresponding lattice size in the physical unit, (1.93 fm) 3 x 3.87 fm). In 
our calculation, the hopping parameter for the s-quark is kept as k s = 0.13710, while the 
three gauge ensembles, K U) d = 0.13760 (Set 1), 0.13800 (Set 2) and 0.13825 (Set 3), are taken 
for u, d-quarks. 
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Table 2 Lattice parameters and hadron masses in unit of [MeV] are listed. 



Lattice parameters 

csw lattice size a [fm] 

L [fm] 

1.83 0.13710 

1.7610 16 3 x 32 0.1209 

1.93 



Nconf 

^ud 

m-Tf 

m K 

m N 

m A 


m s 

Set 1 

700 

0.13760 

875(1) 

916(1) 

1810(2) 

1839(2) 

18466(2) 

1872(2) 

Set 2 

800 

0.13800 

749(1) 

828(1) 

1619(2) 

1675(2) 

1689(2) 

1737(2) 

Set 3 

800 

0.13825 

660(1) 

768(1) 

1482(3) 

1556(3) 

1575(3) 

1640(2) 


Wall source operators which generate positive parity two-baryon states with flavor 
structures h\ and /12 are given by 

^a/3 = [ e abc{QaClbQb)Qca\ h2 [^defiQdC^QDQ fp\ hl , ( 19 ) 

where Q = <z(®) is the quark wall-source. Projection operators for spin-singlet and spin- 

triplet states are given by 

p &°- “ d p «e' = ^4^- (20) 

Quark propagators are calculated for the wall source at to with the Dirichlet boundary 
condition in the temporal direction at t = 16 + to- The wall source is placed at 32 different 
values of to on each gauge configuration, in order to increase the statistics, in addition to 
the average over forward and backward propagations in time. The Af projection of the 
cubic group is taken for the sink operator to obtain the relative S- wave in the BB wave 
function 2 . Numerical computations have been carried out using the KEK supercomputer 
system, Blue Gene/L, and the kaon and jpsi clusters at Fermilab. Hadron masses obtained 
in our calculation are given in Table 2. Thresholds of two-baryons with the strangeness 
S = —2 for each set of gauge configurations are plotted in Fig. 1. 

5. Numerical Results 

We now present our results of coupled channel BB potentials in strangeness S = —2 sector. 
5.1. Time dependence 

We first show how the time-dependent method extended to the coupled channel system 
works in our calculation. For this purpose, we investigate time dependences of the diagonal 
potentials. Fig. 2 shows in 1 <S'o (I = 2) channel (5th line in Table 1) and V n ~ne in 

3 Si (I = 0) channel (2nd line in Table 1) at three values of t — to (= 8, 9,10) with Set 3, which 
has the lightest pion mass in our calculation. Within statistical errors, no significant t — to 
dependence is observed for these single channel potentials with Set 3, showing that t — to = 8 
is large enough to suppress inelastic contributions and that higher order contributions in the 
derivative expansion are negligible. 


2 In this paper, relative D-waves in spin-triplet channels are not explicitly considered but their 
effect is included implicitly in the effective central potentials for spin-triplet channels. 
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Fig. 1 Thresholds of two-baryons with the strangeness S = —2 for each gauge ensemble. 
(Left) A sum of masses for each channel in units of MeV. (Right) A sum of two-baryon 
masses in each channel minus the average of three channels, (m jv + 2m-A + m~ + 2ms)/3 • 




Fig. 2 ££ potential in the 1 S'o (7 = 2) channel (left) and NE potential in the 3 S\ (7 = 0) 

channel (right) as a function of r at t — to = 8(red), 9 (blue) and 10 (green) calculated with 
Set 3. 


In Fig. 3, two diagonal potentials in x 5o (7 = 1) channel (3rd line in Table 1) calculated 
with Set 1 are shown at t — to = 8 ~ 10. Again no significant t — to dependence is observed 
at this quark mass and this is true at other quark masses. Similarly, three diagonal potentials 
in 3 Si (7 = 1) channel (4th line in Table 1) and those in 1 So (7 = 0) channel (1st line in 
Table 1) show no significant t — to dependence at all quark masses, as seen in Fig. 4 for the 
Set 2. 

Since no significant t — to dependences are observed for all diagonal potentials at t — to = 
8, 9,10, we hereafter consider results at t — to = 8, where statistical errors are smallest. 

5.2. Hermiticity 

Hermiticity of the potential matrix is a sufficient condition for the probability conservation, 
though it is not a necessary condition. In this subsection, we investigate the Hermiticity of 
the potential matrix, V a b = V b a , since it is not automatically guaranteed in the definition 
of the coupled channel potential matrix in eq. (5). As in the case of the diagonal parts, we 
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Fig. 3 Diagonal parts of potential matrix in the 1 S'o (I = 1) channel, V n=a n^ (left) and 
F AS AS (right), at t — to = 8(red), 9 (blue) and 10 (green) calculated with Set 1. 
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Fig. 4 (Upper) Diagonal parts of potential matrix in the 3 S\ (/ = 1) channel, V N ^]\ r~ 
(left), D as as (center), and U ss S s (right), at t — to = 8(red), 9 (blue) and 10 (green) cal¬ 
culated with Set 2. (Lower) Same as above but in the x 5o ( I = 0) channel, D AA aa (left), 
V N=L ]\r'E (center), and U ES ss (right). 


confirm that off-diagonal parts of potential matrix show no significant t — to dependence, so 
we take results at t — to = 8 in our analysis. 

We introduce an Hermiticity measured’14_5 = 2 (V a b — V b a )/(V a b + V b a ) to see the relative 
magnitude of the hermiticiy violation of the potential matrix. Fig. 5 presents dU/v--AS in 
the 1 S’o (/ = 1) channel with Set 1 (red), Set 2 (blue) and Set 3 (green). It satisfies the 
Hermiticity well within the statistical errors. 

Fig. 6 shows 5V a _b f° r a,b = in the 3 Si (I = 1) channel. Some violations of 

Hermiticity can be seen in 5V vs-AS and dVAs-ss at r < 0.5 fm region. Those for a,b = 
AA, JVE,££ in 1 S'o (I = 0) are given in Fig. 7. Hermiticity is more or less satisfied within 
the statistical errors. It is our future problem to check whether possible Hermiticity breaking 
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Fig. 5 Hermiticity measure for off-diagonal elements of potential matrix, <5V/vs~ae i n the 
(I = 1) channel on Set 1 (red), Set 2 (blue) and Set 3 (green). 





Fig. 6 Violation of Hermiticity in the channel 3 S\ (I = 1): (left) d'Vvs-AE (center) 
SVns -ee (right) <5Vas-ee on Set 1 (red), Set 2 (blue) and Set 3 (green). 




Set3 — — 

Se'2 ■ 

||ll 




Fig. 7 Violation of Hermiticity in the channel 1 So (1 = 0): (left) <5Vaa-./ve (center) 
<5Vaa-ee (right) ^Vjvs-ee on Set 1 (red), Set 2 (blue) and Set 3 (green). 

for small r in Fig.6 and large r in Fig.7 disappears or not by removing our assumption on 
the wave-function renormalization factor introduced after eq. (1). 

5.3. Potential matrices and their quark mass dependence 

We here separately discuss properties of potentials in three cases, single channel, two channels 
and three channels. 

5.3.1. Single channel. Fig. 8 shows quark mass dependences of the EE potential in 
1 Sq(I = 2) channel (left) and the NE potential in the 3 Si(I = 0) channel (right). We first 
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Fig. 8 The SS potential in the 1 Sq(I = 2) channel (left) and the N E potential in the 
3 Si(I = 0) channel (right). Result from three gauge ensembles, Set 1 (red), Set 2 (blue) and 
Set 3 (green), are shown in one figure. Insets show the enlargement of the same plot. 


notice non-smooth behaviors as a function of r at large r for both cases, which indicate 
that the spatial volume is not sufficiently large. In addition, non-smooth behavior at short 
distance observed in the NE 3 S\(I = 0) channel may be caused by the finite lattice spac¬ 
ing effect. With these systematics, we discuss only qualitative features of potentials in this 
report, and leaves quantitative analysis such as the extraction of scattering phase shifts for 
future studies with larger and finer lattices. 

The SS potential in the 1 Sq(I = 2) channel (left), which belongs to the 27-plet irreducible 
representation in the flavor SU(3), has repulsion at short distance and attraction at long 
distance. Also, the magnitude of these two components increases as the light (ud) quark 
mass decreases, as in the case of the NN potential in the 1 So sector belonging to 27-plet. 
An increase of attraction at long distance, r > 0.8 frn, may be related to the decrease of the 
mass of the pion exchanged between two S’s. 

Similarly, the NE potential in the 3 Si(/ = 0) channel (right) has both repulsion at short 
distance and attraction at long distance. The magnitude of these two components is enhanced 
as the light quark mass decreases. It should be remarked that the repulsion at short distance 
here is weaker than that of the SS potential in the 1 Sq(I = 2) channel. This difference may 
be related to the fact that the Pauli blocking in the quark level for NE in the 3 S'i(/ = 0) 
channel is weaker than SS in the 1 So(I = 2) channel. 

5.3.2. Two channels. A potential matrix in the 1 S'o (I = 1) channel, which has NE and 
AS components, is given in Fig. 9, which shows that diagonal elements of the potential 
matrix in this channel, V n ~ns and F AS A£) are both strongly repulsive and the off-diagonal 
element, W^ae, is comparable to or even larger than the diagonal elements. These features 
have been observed already in the flavor SU(3) symmetric limit [8]. 

5.3.3. Three channels. Fig. 10 shows a potential matrix in the 3 l S'i (/ = 1) channel, which 
has NE, AS and SS components. All diagonal elements of the potential matrix, V n ~ne (left 
upper), F as ae (left middle) and F ss ee (left lower), have an attraction at long distance and 
a repulsive core at short distance. The largest attraction in this channel appears in F ss v;£, 
whose maximum depth is about —60 MeV at around r ~ 0.6 fm. All diagonal potentials 
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Fig. 9 Diagonal (left and center) and off-diagonal (right) elements of the potential matrix 
in the 1 Sq(I = 1) channel. Result from three gauge ensembles, Set 1 (red), Set 2 (blue) and 
Set 3 (green), are shown in one figure. Insets show the enlargement of the same plot. 


have a tendency that magnitudes of both repulsion at short distance and attraction at long 
distance increase as the light quark masses decrease. 

For the off-diagonal elements of potential matrix, F^ae and F AS ee are much smaller 
than F^se- These off-diagonal potentials, F^ae and F ae se, almost vanish at r > 1.2 fm 
and have a small quark mass dependence, while F^se; increases as the light quark masses 
decrease. 

Fig. 11 shows the potential matrix in the 1 5o (/ = 0) channel, where the H dibaryon state 
may appear. All diagonal elements of the potential matrix have a repulsive core at short 
distance, whose strength, however, depends strongly on the state. An attractive pocket, on 
the other hand, appears only in two diagonal elements, F AA aa and V n=l ne, where V Nr 'NE. 
has much deeper attractive pocket than F aa aa does, while F ss se is totally repulsive in 
the whole range of r. 

The off-diagonal element, F AA at~, is smaller than other two, so that the decay rate from 
NS to AA may be relatively suppressed. Diagonal elements of the potential matrix gener¬ 
ated with the configuration Set 2 are most attractive, while off-diagonal potentials with the 
configuration Set 3 are strongest in magnitude for r > 0.5fm. 

5-4- Potential matrix in the SU(3) irreducible representation basis 
We here present potential matrices in the SU(3) irreducible representation basis (SU(3) 
basis in short hereafter) such as 1,8 S , 8 a , 10,10, 27, obtained from the particle basis by using 
Clebsh-Gordan coefficients. This makes us possible to compare the results with those in the 
flavor SU(3) symmetric limit [8]. 

The transformation from particle basis to SU(3) basis is defined as 
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Fig. 10 Diagonal (left three panels) and off-diagonal (right three panels) elements of the 
potential matrix in the 3 S\(I = 1) channel. Result from three gauge ensembles, Set 1 (red), 
Set 2 (blue) and Set 3 (green), are shown in one figure. Insets show the enlargement of the 
same plot. 


for 3 Si(I = 1) channel. 

Fig. 12 shows the potential matrix in the SU(3) basis for 3 S\(I = 1), which is composed 
of 8 a , 10 and 10. While all diagonal elements of the potential matrix have a repulsive core, 
the height of the repulsive core in F 8a 8 a is much lower than other two and its depth of the 
attractive pocket is the deepest among three. On the other hand, F 10 io is strongly repulsive 
and has a quite shallow attractive pocket at all quark masses, though the height of the 
repulsive core and a range of attractive pocket increase as the ud quark mass decreases. As 
far as off-diagonal elements are concerned, they are very small. In particular, F 10 ^ vanishes 
at all quark masses including Set 3, where the SU(3) breaking by the difference between ud 
and s quark masses is maximal in our calculation. Other two off-diagonal elements, F 8 “io 
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Fig. 11 The potential matrix in the 1 Sq(I = 0) channel, as in Fig. 10. 


and F 8 “yq, have small non-zero values at short distance region (r < 0.6 frn), which gradually 
increase as the ud quark masses decrease. 

Fig. 13 shows the potential matrix in the SU(3) basis for 1 5o (/ = 0). As in the case of the 
SU(3) limit [8], the diagonal element for the flavor singlet state, V 1 1, is strongly attractive, 
while is repulsive, at all distances. The absence of repulsive core in V 1 \ is consistent 

with the absence of the quark Pauli blocking effect. A shape of V 27 27 is similar to the 1 Sq 
nuclear force, which also belongs to 27-plet. 

Quark mass dependences of diagonal potentials can be seen clearly in the flavor basis. 
As the light quark mass decreases, the attraction in V x i gradually increases, while both 
repulsive core and attraction at long distance in F 27 27 are enhanced. 

Off-diagonal elements of potential matrix in the SU(3) basis are presented in right three 
panels in Fig. 13, which give effective measures of the flavor SU(3) breaking effects since 
they are absent in the flavor SU(3) symmetric limit. Fig. 13 shows that F 1 8„ (upper) is 


13/17 







































































3000 

2500 

2000 

1500 

1000 

500 

0 

3000 

2500 

2000 

1500 

1000 

500 

0 




0 0.2 0.4 

0.6 0.8 1 
r [fm] 

1.2 1.4 1.6 


. Set3 

^ Set2 

r io Setl 


v 1 


10 


100 
50 
0 
-50 
: -loo 

i 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

i. 






0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

r [fin] 


3000 
2500 
2000 
| 1500 

> 1000 



500 

0 


: -loo 


i 0 0.2 0.4 0.6 0.8 1 

i. 


1.2 1.4 1.6 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 





r [fm] 


r[fm] 


Fig. 12 Potential matrix in the SU(3) basis for 3 Si with S = —2 and 1 = 1. Left three 
are diagonal elements, F 8 “s a (upper), F 10 jq (midle) and F 10 io (lower), while right three 
are off-diagonal ones, F 8 “io (upper), F 8q jq (middle) and (lower). Red, blue and green 

symbols stand for results with Set 1, 2 and 3, respectively. 


small but non-zero while V 1 27 (middle) and F 8s 27 (lower) are consistent with zero except 
for very short distance, r < 0.2 fm, where cutoff effects could be sizeable. These results tell 
us that flavor SU(3) breaking effects in the off-diagonal parts is much smaller than that in 
the diagonal part within the quark masses adopted in this paper. The 27-plet state is almost 
uncoupled even if ud quark mass is different from the s quark mass. 

For 1 5 q(I = 1), the potential matrix in the SU(3) basis is obtained by 
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Fig. 13 Potential matrix in the SU(3) basis for 1 Sq with S = —2 and 1 = 0. Left three 
are diagonal elements, V 1 \ (upper), F 8s 8 s (middle) and F 27 27 (lower), while right three 
are off-diagonal ones, F 1 ^ (upper), V 1 27 (middle) and F 83 27 (lower). Red, blue and green 
symbols stand for results with Set 1, 2 and 3, respectively. 


Fig. 14 shows diagonal and off-diagonal parts of potential matrix in the SU(3) basis for 
1 Sq(I = 1) channel. We find that diagonal elements, F 8 ^ and F 27 27 , have similar behaviors 
to ones obtained from 1 Sq(I = 0), and the transition potential between 8 s -plet and 27-plet 
is quite small for all Sets. 

To see the effects of SU(3) breakings, it is interesting to compare the potentials in the 
SU(3) basis extracted from two different channels, 1 Sq(I = 0) and 1 So(I = 1). In Fig. 15, we 
show the measure of SU(3) breaking defined as AF ^ = 2 (Vf =0 — Vf =1 )/(Vf =0 + Vf =1 ). We 
find that, for AI/( 8s ), there are no significant deviations from zero within statistical errors 
for all Sets. A similar behavior to AF( 8 d case can be seen again for <5F^ 27 ^ case except for 
a singular behavior at r ~ 0.4 fm where the 27-plet potentials themselves almost vanish. 
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Fig. 14 Diagonal (left and center) and off-diagonal (right) elements of the potential matrix 
in the SU(3) basis for the 1 Sq(I = 1) channel. Left and center panels are diagonal elements, 
F 8 s 8 s (Left), F 27 27 (center), while right panel is off-diagonal ones, V 8s 27 - Red, blue and 
green symbols stand for results with Set 1, 2 and 3, respectively. 




Fig. 15 Measures of SU(3) breaking effects: (left) and SV^ (right) as a function 

of r. Red, blue and green symbols stand for results with Set 1, 2 and 3, respectively. 


6. Summary and conclusions 

In this paper, we have investigated the S = —2 BB potentials from 2 + 1 flavor lattice QCD 
by using the HAL QCD method extended to coupled channel systems in Ref. [9]. Combining 
the coupled channel formalism with the time-dependent Schrodinger equation [18], we could 
extract potential matrices for the first time without the ground state saturation and without 
the diagonalization of the source operators. 

By considering two baryon systems with S = — 2, AA, NS, EE and AS which are mutually 
coupled, we successfully extracted potential matrices. They are approximately hermitician 
within the statistical errors, which is not guaranteed from the definition. A small violation 
of Hermiticity may be removed at least partly by the proper treatment of renormalization 
factors, which is left for future studies. 

We discussed properties of potential matrices for all S = —2 two baryon systems. We found 
that all diagonal elements of the potential matrix have a repulsive core, while their heights 
largely depend on their flavor structure. Our previous works show that decreasing ud quark 
mass leads to the enhancement of the short ranged repulsion and the long-ranged attraction. 
Although such quark mass dependence was seen clearly in single channel case, it becomes 
less pronounced in two and three channel cases. 
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The potentials in the SU(3) basis are also investigated, where we could see clear quark 
mass dependence. We found a strongly attractive potential for V 1 \ whose strength increases 
as the quark mass decreases. The off-diagonal potentials in the SU(3) basis is a proper 
measure of the SU(3) breaking. In the 3 Si(/ = 1) channel, except for W 10 io, only a small 
transition potential between irreducible representations could be seen at short distances. In 
the 1 S' 0 (I = 0 ) channel, we found a clear mixture of the flavor singlet state and the octet 
state. The other off-diagonal potentials have only a small magnitude at short distances. Note, 
however, that the SU(3) breaking introduced in this paper ( m n /rriK = 0.96,0.90,0.86) is still 
small compared to the realistic magnitude of the breaking {m^/mK = 0.27). Nevertheless, 
the present paper provides a first theoretical and numerical step toward the realistic BB 
potential matrix at the physical quark masses. 
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